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Stabilization of photon-number states via single-photon corrections 
a first convergence analysis under an ideal set-up 


H. B. Silveira P. S. Pereira da Silva P Rouchon 


Abstract —This paper presents a first mathematical conver¬ 
gence analysis of a Fock states feedback stabilization scheme via 
single-photon corrections. This measurement-based feedback 
has been developed and experimentally tested in 2012 by 
the cavity quantum electrodynamics group of Serge Haroche 
and Jean-Michel Raimond. Here, we consider the infinite¬ 
dimensional Markov model corresponding to the ideal set¬ 
up where detection errors and feedback delays have been 
disregarded. In this ideal context, we show that any goal Fock 
state can be stabilized by a Lyapunov-based feedback for any 
initial quantum state belonging to the dense subset of finite 
rank density operators with support in a finite photon-number 
sub-space. Closed-loop simulations illustrate the performance 
of the feedback law. 

I. INTRODUCTION 

In [8], a photon-number states (Fock state) feedback stabi¬ 
lization scheme via single-photon corrections was described 
and experimentally tested. Such control problem is relevant 
for quantum information applications [6], [4]. The quantum 
state p corresponds to the density operator of a microwave 
field stored inside a super-conducting cavity and described as 
a quantum harmonic oscillator. At each sample step fc G N, 
a probe atom is launched inside the cavity. The measurement 
outcome yk detected by a sensor is the energy-state of this 
probe atom after its interaction with the microwave field. 
Each probe atom is considered as a two-level system: either 
it is detected in the lowest energy state \g), or the highest 
energy state |e). Consequently, the measurement outcomes 
corresponds to a discrete-valued output yk with only two 
distinct possibilities: g or e. Similarly, the control inputs 
Uk are also discrete-valued with 3 distinct possibilities: 
— 1,0,4-1. The open-loop value Uk = 0 corresponds to 
a dispersive atom/field interaction: it achieves in fact a 
Quantum Non-Demolition measurement of Fock states [2]. 
The two other values = ±1 correspond to resonant 
atom/field interactions where the probe atom and the field 
exchange energy quanta: these values achieve single-photon 
corrections. 
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Although the feedback law proposed and implemented 
in [8] considered imperfect detections on yk and delays 
in the control, here we focus on an ideal-set up, that is, 
detection errors and control delays have been disregarded. 
Theorem|2] shows that, by adding an arbitrarily small term to 
the Lyapunov function used in [8], one ensures almost sure 
global stabilization of any goal Fock state for the closed- 
loop ideal set-up. This is achieved by relying on an infinite¬ 
dimensional Markov model of the ideal set-up that takes into 
account the back-action of the measurement outcome yk on 
the quantum state pk+i- 

Loosely speaking, in [8], the control value Uk at each 
sampling step k was chosen so as to minimize the con¬ 
ditional expectation of the Lyapunov function V{pk) = 
Ti {d{N)pk), where TV is the photon-number operator, 
d{n) = {n — n)^ and p = |fT)(n| is the goal Fock state. 
However, in closed-loop, the difference between such V 
and its conditional expectation is not strictly positive: such 
V does not become a strict Lyapunov function in closed- 
loop and additional arguments have to be considered to 
prove convergence. These additional arguments are related 
to Lasalle invariance. They are well established in a smooth 
context where the control it is a smooth function of the state 
p. This cannot be the case here since it is a discrete-valued 
control. In order to overcome such technical difficulties, we 
propose, similarly to [1], to add the arbitrarily small term 
l^(Pfc), where e > 0. This slightly 
modified control-Lyapunov function becomes then a strict- 
Lyapunov function in closed-loop that simplifies notably the 
convergence analysis. Moreover, the developed convergence 
analysis is done in the infinite-dimensional setting in the 
following sense: we show that, for any initial density operator 
po with a finite photon-number support (po|n) = 0 for n 
large enough), the closed-loop trajectory k i-A pk remains 
also with a finite photon-number support with a uniform 
bound on the maximum photon-number. This almost finite¬ 
dimensional behavior simplifies the convergence analysis 
despite the fact that such condition on po is met on a dense 
subset of density operators (Hilbert-Schmidt topology on the 
Banach space of Hilbert-Schmidt self-adjoint operators). 

The paper is organized as follows. Section [H] presents 
the ideal Markov model of the experimental set-up of the 
controlled microwave super-conducting cavity reported in [8] 
and precisely formulates the Fock state stabilization problem 
here treated (see Definition [T]). Section |III] establishes the 
proposed solution to the control problem in two distinct parts. 
Firstly, Section IIII-AI considers the case where the initial 
condition po is a diagonal density operator (see Theorem[T]). 




Only the main ideas of the convergence proof are outlined. 
The technical details are given in Section |Vl Afterwards, in 
Section IIII-BI the main result of the paper is presented: the 
general solution is obtained from Theorem[T]for po belonging 
to a dense subset (see Theorem |2|. The simulation results 
are exhibited in Section |IV] The proof of some intermediate 
results and computations required in Sections |III] and |V] 
are presented in Appendices IB14 g 1 Finally, the concluding 
remarks are given in Section |Vll 


II. IDEAL MARKOV MODEL 


Denote by Ti, the separable complex Hilbert space L 2 (C) 
with orthonormal basis {|n), n G N} of Lock states (photon- 
number). Hence, U = (V'o, V’l: • ■ •) G 

(2(C)}. Let D be the set of all density operators on T-L, that 
is, the set of trace-class, self-adjoint, non-negative operators 
on Ti. with unit trace. The sample step, corresponding to 
a sampling period around 100/is, is indexed by fc S N = 
{0,1,2,...}, Uk G { — 1,0,1} is the control, G D the 
quantum state and £ {g, e} the measurement outcome. 
The ideal Markov model of the controlled microwave super¬ 
conducting cavity used in [8] is given by: 


Pk+l — 


g _ M g{uk)pkM\{uk) 
^'=+1 “ l>[{Mg{uk)pkMl{uk)) 
e _ Me{uk)pkMl{uk) 
^'"+1 “ Tv{M,(uk)pkMl{uk)) 


when yk= g, 
when yk = e, 

( 1 ) 


where the measurements outcomes yk = g and yk = e 
occur with probabilities pg^k — Tr (^Mg{uk)pkM^g{uk)^ 

and = Tr ^Me(ufc)pfcM|;(Mfe)^ = 1 - Pg,fc, respec¬ 
tively, Uk = 0 corresponds to a dispersive interaction of 
the launched atom with the cavity field (Quantum Non- 
Demolition measurement of photons) 


Mg(0) = COS , Me(0) = sin , (2) 


when Uk = +1 the atom enters the cavity in the state |e) 
with a resonant interaction with the cavity field 


Mg( + 1) -- 

when Uk = 



Me( + l)=COS , 

(3) 


1 it enters in \g) with a resonant interaction 


Mgi-1) = cos , Me(-l) = a \ 

(4) 

and (j)o,4>R,0o G R are adjustable control parameters. Lor 
each u G { — 1,0,1}, Mg{u) and MeM are (linear) 
operators on "H defined in the obvious wa}|3 according to 
the definitions in Appendix lAl They are indeed well-defined 
operators on H, despite the fact that a and are unbounded 


*As usual in quantum physics, it is here assumed that the measurement 
outcome yp = V cannot occur when Tr (^My{uk)pkMy(uk)^ = 0, for 
y = 9,e- 

^For instance, Mg{+l)\n) = (^sm{^\fN)/\/~N^ yjn -|- l|n -|- 1) = 

sin(^ ^n+ l)|n + 1). In order for the definition of Me(—1) to be 
consistent, it is assumed sin(0)/0 = 1. 


operators. It is clear that Mg{u), Me{u) are bounded oper¬ 
ators on V. with M^g{u)Mg{u) + MI{u)M e{u) = I (iden¬ 
tity operator), = AT^(-l-l) = asin(^-\/iV)/-\/iV, 

and Afg(—1), A4’g(0), Afe(O), Afe(+1) are self-adjoint. It 
is easy to see that if the initial condition po is a density 
operator then, for all realizations of the ideal Markov pro¬ 
cess O-©, Pk is a density operator for fc G N. 

Notice that p = |n)(rl| is a steady state of the Markov 
process ©-© with Uk = 0, where n G N is arbitrary. The 
control problem here treated is given as follows: 

Definition 1: Lor the ideal Markov process ©-©, the 
control problem is to find a feedback law Uk = f{pk) such 
that, given an initial condition po and n G N, the closed-loop 
trajectory pk converges almost surely towards the goal Lock 
state p = |fi)(fl| as k ^ 00 . 

The almost sure convergence above is with respect to the 
probabilities amplitudes Pn{p) = Tr(|n)(n|p) = (n|p|n) of 
p, that is, limfc_,.oo Pn{pk) = Pn{p) for each n G N. In other 
words, limfc^oo Pn{pk) = 1 and limfc_>oo Pn{pk) = 0 when 
n ^n. The solution proposed in this paper for the control 
problem above is developed in the next section. 

HI. STABILIZATION OL LOCK STATES 

Given any operator A: TL —>■ H, let Amn = {m\A\n) for 
TO, n G N. Hence, is the n-th diagonal element of A, 

while Amn with m n correspond to its “off-diagonal” 
elements. One says that the operator A is diagonal when 
Amn — 0 for all to, n G N with m ^ n. One shall begin 
by solving the control problem given in Definition [T] in the 
particular case where the initial condition po is diagonal (see 
Theorem [1] in Section IIII-Al i. Afterwards, in Section IIII-BI 
the solution to the general non-commutative case is presented 
(see Theorem |2]i: its solution relies essentially on the diago¬ 
nal case. 

A. Diagonal case 

Lor each n* G N, defin^l 

Dn* = {p G D I p is diagonal and p\n) = 0,Vn > n*} . 

Consider the set Note that C 

Dn*+i, and that each element p of D, is “finite dimensional” 
in the following sense: p G D is in if and only if 
p = P'nn\'kt){n\, and p G may be considered as 

an operator from TL to the finite-dimensional space "Hn* = 
span{|0),..., |n*)}, or as a density matrix on Tin*- One 
defines the functions rimin' —>■ N, rimax- —>■ N and 

niength'- D* N respectively by: 

• rimin{p) is the smallest n G N such that p\n) ^ 0; 

• rimaxip) i^ ths greatest n G N such that p\n) ^ 0; 

• rii^ngthi^p) — ttmaxip') rimini^p)- 

It is clear that, given p G Dt,, one has p G if and only 
if rimaxip) ^ ri*. The next result exhibits the properties of 
the state pk of ©-© with respect to these functions. 

^Note that if p = |n){n| for some n G N, then p G Dn- 

















Proposition 1: For every realization of the ideal Markov 
process ([T])-® with initial condition po G 13*, one has that 
Pk G for all fc G N with: 

• If tik = 0 or Ilk — 1, then rij^axipk+l) ^ ^max(PA:) 

and Tilengthi^Pk+l) ^ tliengik^Pk)) 

• If Uk — “fl, then vT-max(pfc+i) ^ tijYiaxiPk) 1 and 
ttlengthiPk+l) — ’tllength(yPk'} ■ 

Proof: See Appendix IbI ■ 

Take a goal photon-number n G N. As in [1], consider the 
following Lyapunov function I^: Z3, —K defined as 

Ve{p) = Tr {d{N)p) - e ^ pL: for p G £>*, (5) 

raSN 

where e > 0 is a real number and d{n) = (n —n)^ as defined 
in [8]. The feedback law u: Z3* —{—1,0,1} is given by 

u = f{p) - Argmin E [V^{pk+i) \ pk= p,Uk = v]. (6) 

Note that for each p G -D* and n* > nmax{p )7 d{N)p in 
© is a well-defined self-adjoint, non-negative, trace-class 
operator on H, by considering d{N) as an operator on TLn* 
and p as an operator from TL to TLn*- Indeed, d{N)p = 
— n)^|n)(n|. Thus, © is well-defined. More¬ 
over, since "Hn* is invariant under p G for n* > rimaxip), 
it is clear that Ty {d{N)p) = Tr^^, (d(7V)p), where on the 
right-hand side one considers p as an operator on the finite¬ 
dimensional space Hn* and the trace is taken over TLn* ■ 

We have the following convergence result when po G H*: 
Theorem 1: Let n G N and e > 0. In ©-©, assume that 
fo/Tt and (^o/tt)^ are irrational numbers, and take (j>R = 
-K12 — nfo. Consider the closed-loop Markov process ©- 
© with Uk = f{pk), where the feedback law / is as in ©. 
Then, given any initial condition po G H*, one has that pk 
converges almost surely towards p = |n)(n| as A: —>■ oo. 

Its proof is decomposed into two steps: 

First Step. Choose Tt G and c ^ 0. Let tiq — 
ro = n-mm(po)- Then, there exists an integer mg >no + 
ro -f n -f 1 (depending on no, ro, n and e) such that, for all 
closed-loop realizations pk, one has pk G Dmo for fc G N. 
Second Step. Choose irrational numbers fg/Tt and {Og/Ttf^ 
in and take fn = it 12 — nfg. In DmQ, 14 is a strict 

super-martingale: for all density operators p in Dmo 7 0 *^^ 

E[V',(pfe+i) I Pk = p,uk = /(p)]-l4(p) = -Qk(p,/(p)), 

where Qv^ (p, /(p)) > 0, and Qk (p, /(p)) = 0 if and only if 
p = 'p. The almost sure convergence follows then from usual 
results on strict super-martingales for Markov processes with 
compact state spaces. 

The complete proof of the two steps above is presented in 
Section lY] The general case where the initial condition pg 
is not necessarily diagonal is treated in the next subsection. 

B. General case 

Consider, for each n* G N, 

®n* = (p G D I p\n) = 0,Vn > n*} C ID„*+i, 


and let D* = Un*GN®"* ^ clear that p G D is in 

D„. if and only if p = Consequently, 

D* is a dense subset of D when D is endowed with the 
subspace topology induced from the Hilbert-Schmidt norm. 
Indeed, let JT 2 be the complex Banach space of all Hilbert- 
Schmidt operators on "H with the Hilbert-Schmidt norm 
\\B \\2 = for B GJ 2 [7], [3]. Since 

D C J72 and p G Dn. has the form p = Em.n=o Pmn\m){n\, 
the density property of D* in D is clear. 

One has that p G Dn* may be considered as an operator 
from TL to the finite-dimensional space TLn *, or as a density 
matrix on TLn*- Hence, d{N)p is a well-defined trace-class 
operator on "H, by considering d{N) as an operator on 
Tin* and p G as an operator from TL to "Hn*- Indeed, 
d{N)p = Em,n=0 (to —n)^|m)(n|, and it is trace-class 
because its range is finite-dimensional [7], [3]. Consequently, 
the Lyapunov function 14 in ©, the feedback in © and 
timax can be extended to IJ,. 

Define the map A: D* —>■ 13* C D* as Ap = 
Px^„\n){n\. Note that A extracts the diagonal of 
p G D*. It is easy to see that Umaxi^p) = timaxip) and 
(Ap)„„ = pnn 7 p G D*. Moreover, Ap = p when p G D^. 
Other properties of the map A are given in the next result: 

Proposition 2: Let p G D*, u G {—1, 0,1}, y = g,e. Take 
a = Tr {My{u)pMl{u)). Then: 

• Tr (Ap) = Tr (AAp), for every diagonal bounded 
operator A: TL ^ TL; 

• 14(p) = I4(Ap), for e > 0; 

• a~^My{u)pM^y{u) belongs to ID* with 

A{a~^ M y{u) pM^y{u)) = a~^ M y{u){Ap)M^y{u); 

• [My{u){Ap)Ml{u)]^^ = [My{u)pMl{u)]^^, for 

all n G N. In particular, a = Tr (^My(u){Ap)M^y{u^. 
Proof: See Appendix iGl ■ 

Now, let e > 0 and p = |n)(n|, where n G N. Assume that 
Po G D*. Let pfc, fc G N, be the corresponding closed-loop 
trajectory for a fixed realization of ©-© with feedback 
Uk = /(pfc), where / is as in ©. It is immediate from the 
proposition above that: 

• Pk G D*, for fc G N; 

• Apk G 13*, A; G N, is the corresponding closed-loop 
trajectory of ©-© for the initial condition Apo, the 
same realization (and with the same transition proba¬ 
bilities pe,fc and Pg,fc), as well as the same feedback 
Uk = f{pk) = /(Apfc); 

• Tr (|n)(n|pfe) = Tr (|n)(n| Ap^), for any n G N. 

From these arguments. Theorem [T] and the fact that Ap = 

p, one immediately obtains the following generic solution 
to the control problem, that is, when the initial condition pg 
belongs to the dense subset D* of D: 

Theorem 2: Let n G N and e > 0. In ©-©, assume that 
fg/it and {Og/Ttf^ are irrational numbers, and take (j>B. = 
1:12 — rifg. Consider the closed-loop Markov process ©- 
© with Uk = /(pfc), where the feedback law / is as in ©. 
Then, given any initial condition pg G D*, one has that pk 
converges almost surely towards p = |n)(n| as A: —>■ oo. 


IV. SIMULATION RESULTS 


Goal Fock state = 10 


This section presents the closed-loop simulation results 
concerning the application of Theorem |2] above to the 
ideal Markov process ([T]i-(l4]i. The quantum experimental 
results exhibited in [8] used the following control parameter 
values in (|2|i-(|4|i: (/iq/tt = 0.252 and ^q/tt Ri 2!\Jn + 1. 
However, according to the assumptions in Theorem |2l (/iq/tt 
and should be irrational numbers. Hence, here one 

chooses (/)o/3.14 = 0.252 and 0o/3.14 = 2l^n + 1. One 
takes po = Sn=o \'n){n\/\& £ D* as the initial condition, 
n = 10 for the goal Fock state p = |n)(n|, and e = 10^ as 
the gain for the feedback Uk = f{pk) in dU-®. Figure [T| 
exhibits the simulation results for one closed-loop realization 
with such choices and a final sample step of 120. It shows: 
the dynamics of the populations of pk (top), the controls 
Uk (middle) and the simulated outcomes yk (bottom). The 
population^ of pk correspond to the following observables: 
^1 = J2n=o\n){n\ (n < n), A 2 = \n){n\ {n = n), 
A 3 = (t^ ^ Therefore, one sees from 

the dynamics of the populations that pk converges to p as 
k 00 , which is in accordance with Theorem |2] Note that 
{n\pk\n) ~ 1 and Uk = 0 for all k > 45. 

Recall that Theorem |2] assumes that e > 0. In order 
to further analyze the performance of the Lyapunov-based 
feedback law here proposed, we now make a comparison 
with the one used experimentally in [8], which corresponds 
to take e = 0 in 0, i.e. to disregard the term —e Pnn- 

Figure |2] presents the simulation results for one closed-loop 
realization of such case. The control parameters, po and n = 
10 are the same as above. Note that (n|pfc|ri) Ri 1 and Uk = 0 
for all k > 78. In order to make a comparison in terms of 
the speed of convergence, define the settling time kg to l^e 
the smallest k € N such that (njpkln) > 0.9 for all k > k. 
One has kg = 45 for the case e = 10^ above, and kg = 78 
for e = 0. Therefore, in the two realizations here simulated, 
the choice of e = 10^ reduced the settling time kg by nearly 
42% with respect to e = 0. This behavior is typical on an 
average basis, thereby justifying the term —iJ2nenPnn 
0. Table H shows the average value kg and the standard 
deviation a of kg for e £ {0,0.1,1,10,10^, 10^, 10"^, 10^}, 
where a total of 5000 realizations were simulated for each 
e. Notice that when e is relatively large or relatively small 
in comparison to e = 10^, the average settling time kg 
deteriorated. Furthermore, although for e = 10^ one has that 
kg increased by nearly 22% in comparison to e = 10^, the 
standard deviation a decreased by nearly 62%. Computer 
simulations have suggested that a choice of e > 0 which 
may perhaps significantly improve kg generally depends on 
the initial condition po and on the goal Fock state p = |n)(n|, 
and it has to be determined heuristically. 

V. PROOF OF THEOREM [U (DIAGONAL CASE) 

Proof of the First Step: 

Let e > 0. Define V: —>■ R and W: —>■ R as 

U(p)=Tr(d(Ar)p), W{p)=-J2pL, (7) 

n^N 
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Fig. 1. Simulation of one closed-loop realization with gain e = 10^: 
convergence of towards p (top), controls (middle), and outcomes 
(bottom). Notice that {n\pk\n) ^ 1 and = 0 for all k > 45. 
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Fig. 2. Simulation of one closed-loop realization with gain e = 0: 
convergence of pk towards p (top), controls nj. (middle), and outcomes 
yk (bottom). Notice that {n:|pfc|n} Ri 1 and Uk = 0 for all k > 78. 


TABLE I 

Average settling time kg and standard deviation o- as a 
FUNCTION OF THE GAINS e, CONSIDERING 5000 REALIZATIONS 


e = 0 

kg = 79.94 
o- = 164.97 

e = 0.1 
kg = 79.95 
o- = 166.61 

e = 1 

kg = 81.24 
cr = 174.29 

e = 10 
kg = 71.33 
cr = 150.95 

e = 10^ 

kg = 60.41 
o- = 119.39 

e = lO'^' 

kg = 44.18 
a = 44.12 

e = lO"" 

kg = 47.05 
cr = 37.37 

e = lO'^ 

kg = 53.77 
cr = 16.84 


respectively. Note that Vg = V + eW. Define: 

• Qwip,u) = lU(p) - E [W{pk+i) \ Pk= P,Uk = u], 
. Qv{p, w) = u(p) - E [V(pk+i) I Pk = P,Uk = u], 

. QvSp^u) = Ve{p) - E [v;(pfc+i) I pk = p,uk = u], 

































































































































for p € Di, and u G { — 1, 0,1}. The proof of Theorem [T] is 
a straightforward consequence of the next proposition: 

Proposition 3: Let e > 0 and nQ,ro,n G N. There exists 
an integer toq > no + ro+n + 1 (depending on e,no, tq, n) 
such that, for each p £ with niength{p) < no, if 
rimaxip) = mo, then 


> max{Qy^(p,0),(5y^(p,+l)}. 

In fact, given po G let uq = niength{po) and 

^0 = riminipo)- Note that Umaxipo) = no + ro < mo- By 
Proposition!!] pk G I?* with niength{pk) < no, for all fc G N. 
Since u = f{p) maximizes Qv^{p,f{p)), Proposition [3 
implies that when Umaxipk) = mo for some k £ N, then 
the input Uk will be always be equal to —1, and hence 
Proposition [T] ensures that nmax{pk+i) < n^axipk) = mo- 
Therefore, rimaxipk) < mo, fc G N, showing the First Step. 

The following two lemmas are instrumental for showing 
Proposition [3] Their proofs are given in Appendix |D] and 
Appendix respectively. 

Lemma 1: Given an arbitrary nonzero Oo G K, fix any 
a G M such that 0 < a < 1/2. For all nonzero No, N G N, 
there exists an integer N > N big enough such that, 

0 < 1/2 — a < sin^ (^V^) ^ 1/2 + a, 

for n = N, N + 1, ■ ■ ■, N + Ng — 1. 

Lemma 2: Let p £ D,. Then: 

• \Qw{p,u)\<^, for each u G {—1, 0,1}; 

• Qvip, 0) = 0; 

• Qv{p,+1) = -^p„„[2(n-n) + 1] sin^(^\/n + l); 

kGN 

• Qvip,-^) = '^Pnn [2(n-ri) - 1] sin^ 

The proof of Proposition [3 is shown in the sequel. 
Proof: Let e > 0 and no,ro,n £ N. One has to show that 
there exists mo > ng + tq + n + 1 such that, if p G D^, with 
tiiengthip) < tio, then M = — 1 always maximizes QvSPi'^) 
whenever rimaxip) = mg- From Lemma |2] and the fact that 
Qv, = Qv + tQw, to complete the proof it suffices to show 
that: 

• If p G D:,, is such that ni^rigthi^P^ tig and Tlmax(,p} ^ 
ng + n, then Qv{p,+f) < 0; 

• There exists mg > ng + rg + n + 1 such that 
Qv{p,—f) > 2e, whenever p G H* is such that 
niengthip) < ng and '^max (p) = mg. 

Note that 

'^max (p) 

Qv{p,+l) = - ^ p„„[2(n-n) + l]sin^ (^v/rT+T) , 

n=nmi-n{p) 


for any p £ D,,. Thus, if riiengthip) < ng and Umaxip) > 
n+ng, then Uminip) > n, and hence the first claim is shown. 
Now, fix 0 < a < 1/2 and le0 N > ^ + 2ri + 1 


^l/2-a 

Applying Lemma [T] for Ng = ng + rg + l and such choice of 
N, one gets At > At in which 0 < 1/2 — a < sin^ 


^As is an integer, it follows that N > n 1. 


forn = A^, A^ + 1,..., A^ + no + ro- Take mo = A^ + no + ro- 
Let p £ Dt with niengthip) < tig and nmaxip)_= mg. Note 
that mg > ng + rg + n + 1 and Uminip) > N rg. From 
Lemma |2] and the inequality above for 1/2 — a, one obtains 

mo 

Qv{p,-f)= ^ p„„[2(n - n) - 1] sin^ 

n — rimin (p) 

mp 

> ^ p„„[2(n - n) - l](l/2 - a) 

n—njnin{p) 

mo 

> ^ p„„[2(7V-n)-l](l/2-a) 

n—rijninip) 

mo 

= [2(7V-n)-l](l/2-a) ^ p„„. 

n — rimin (p) 


Using the fact that 


2 e 


l/2-i 


+ 2n + 1 


mo 

n=nmin{p) 


= 1 and N > 


, one shows the second claim, thereby 


completing the proof of Proposition |3| 


Proof of the Second Step: 

Let e > 0. Recall that, by definition, Qv^ — Qv + ^Qw- 
Using the same notation of the First Step, the central idea 
of the proof is to show that, given p G Dmo, one has that 
QvSPj fip)) > 0’ ^nd that QvAP^fip)) = 0 if and only if 
p = p. The following lemma is instrumental for the proof of 
such property. Its proof is presented in Appendix 10 

Lemma 3: Assume that fg/Ti is an irrational number in 
@, and take fa = tt 12 — nfg, where n G N. Let p G H*. 
Then: 

• QwiP^O) > 0, and Qwip, 0 ) = 0 if and only if p = 
|n)(n| for some n G N; 

• Qwip,+l) = Qwip,-^) = 0 whenever p = \n){n\ 
for some n G N. 

One has that mg > n, and (0g/7r)^ is an irrational number 
by assumption. Recall that sin^ (cc) = 0 if and only if cc = ^^t, 
where f is an integer. First we show that Qve(Pj/(p)) = 
0. By Lemma 121 Qvfp,+f) = — siiA(^-\/^ + 1) < 0; 
(5y(p, —1) = — sin^(^v^) < 0 when n > 0, and 
Qvfp,—f) = 0 when n = 0; and Qy(p, 0) = 0. As 
Qwip,u) = 0, for u £ {—1,0,1}, and u = f{p) maximizes 
(5v^(p, u), one has that QvSP^fip)) = 0- 

Now, let p G Dmo C Since u = /(p) maximizes 
(p, u), it follows that 

QvAP^fip)) > Qv(PjO) + eQwip,0) = eQtv(p,0) > 0. 


Suppose QvAPj fip)) = 0- Hence, Qwip,0) = 0, and so 
p = |ri)(n| for some n £ (0,1,..., mg}. It suffices to show 
that QvApJip))\\n){n\ > 0 for n G (0,1,....mg} with 
n Assume that n > n. It is clear that u = fi\n){n\) = 
-land QvApJ{p))\\n){n\ = [2(n-n)-l] sin^(|^v^) > 0. 
Assume now that n < n. Then, u = /(|n)(n|) = +1 and 
QvSpJip))\\n)(n\ = -[2(n-n) + 1] sin^(|iV« + 1) > 0. 
This completes the proof of the referred property. 

The remaining part of the proof of the Second Step is 
a straightforward consequence of the standard stochastic 
convergence result below: 












Theorem 3: [5, Theorem 1, p. 195] Let 17 be a probability 
space and let bL be a measurable space. Consider that 
Xk'. 17 —?> W, k G N, is a Markov chain with respect to 
the natural filtration. Let ® ^ —>■ M be 

measurable non-negative functions with V{Xk) integrable 
for all A: G N. If E [V{Xk+i) \ Xk] - V{Xk) = -Q{Xk), 
for k G N, then limfe_>oo QiXk) = 0 almost surely. 

Indeed, let J'l be the complex Banach space of all trace- 
class operators on H with the trace norm || • ||i, that is, 
||S||i = Tr(|S|), where \B\ = y/WB, for B G Ji. Recall 
that ||i3|| < ||i?||i ancd | Tr (AS) | < ||A||||i?||i, for every 
B G Ji and each bounded operator A: TL ^ TL, where || • || 
is the usual operator norm (sup norm of bounded operators) 
[7], [3]. Consider the subspace topology on Dmo with respect 
to Ji. One has that the closed-loop trajectory pk, k G N, 
is a Markov chain with phase space (with respect to 
the natural filtration and the Borel algebra on Dmo)- R is 
clear that 77^0 is compact, and that Qg and 14 — are 
non-negative and continuous on Dmo, ^ii e > 0, where 
tte = minpg£)^^ Re(p)- The theorem above implies that pk 
converges almost surely towards p as /c —>■ oo (with respect 
to the trace norm). This completes the proof of Theorem [T] 

VI. CONCLUDING REMARKS 

This paper provided a convergence analysis of Fock states 
stabilization via single-photon corrections under an ideal set¬ 
up, that is, assuming perfect measurement detection and no 
control delays. In terms of convergence speed, the simulation 
results here presented have justified the inclusion of the 
term —ill Ih® Lyapunov-based feedback law 
(ID-®. It is straightforward to verify that the convergence 
analysis developed in this paper remains valid for; (i) any 
other function d{n) in Q satisfying d(n) = 0, d{n) is 
increasing for n > n and d{n) is decreasing for n < n; 
and (ii) e > 0 dependent on n, that is, to take the term 
— EraGN ^^Pnn- However, it is an open problem how to 
choose the function d{n) and the gains e„ > 0 so as to 
achieve the best convergence speed. 

Finally, the feedback law used in [8], which corresponds 
to e = 0, was tailored for an experimental set-up with mea¬ 
surement imperfections and control delays. The convergence 
analysis of such realistic situation will be investigated in the 
future. 
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Appendix 

A. Basic properties of the operators N, a and 

Fix n* G N and let Hn* = span{|0),..., |n*)}. Con¬ 
sider the (linear) operators TV; TLn* T~Ln*, a: TLn- 
Tin*-! C Tin*, Cl^- T~Ln* 'Hn*+l defined respectively as 
TV|n) = n \n), a|0) = 0, a\n) = y/n \n — 1) for n > 1, 

^One also recalls that if A is a bounded operator on B. and B S jTi, 
then AB, BA e Ji with Tr (AB) = Tr (BA). 


a^n) = s/n + 1 |n -I- 1). Note that these operators cannot 
be extended to TL. Let /; N —?► M be a function. Define 
the operator /(TV); Hn* -A- Hn* by /(TV)|n) = f(n)\n), 
for each n = 0, ...,n*. It is clear that /(TV) can be 
extented to TL whenever / is a bounded function. Given 
/; N —?> M and an integer m, one defines g; N —?> K as; 
g{n) = f(n + m), when n + m > 0; and g(n) — 0, 
when n + m < 0. One abuses notation letting /(TV + m) 
stand for g{N). Given two functions /, p; N —>■ K., it 
is clear that f(N)g(N) = g(N)f(N) = (fg)(N) and 
(/ -f 5 )(TV) = /(AT) -f g{N). Furthermore; aa^ = N + I, 
a^a = TV, a/(TV) = /(TV -f l)a, a^f(N) = /(TV - l)at. 


B. Proof of Proposition [7] 

Fix any p G T7* and let n G N. In particular, p\n) = 
PnnW)- It then follows from that; 

TVT,(0)pTVTt(0)|n) =p„„cos2 |n), 

Me{0)pM\mn) = pnnsin^ \ri), 

Mf,{Tl)pM\(+l)\n) = Pnn COS^ {^Vn+ l) |n), 

Mg{-l)pMl{-l)\n) = Pnn cos^ (tV^) I^)> 
Me{-l)pM\{-l)\n) = Pn+I,n+1 sin^ {^sfnTT) \n). 
Therefore; 

rimaxip) 

MMpMI(0)= Pnncos" (^2^) |n}(n|, (8) 

rimax (p) 

MA0)pMl(0)^ PnnSin" |n)(n|, (9) 

n=ny^.j.^ (p) 

Timax (p)"i"l 

Mg{ + l)pMl{ + l) = ^2 Pn-l.n-1 Sm^ (10) 

”'=''T'min(p) + l 
rimaxip) 

Me(+l)pMl{+l) = ^2 Pnn cos^ (^vTT+T) |n}(n|, (11) 

n = n^ 271 (p) 

n-maxip) 

Mg{-l)pMl{-l) ^ ^ p„„cos^ (^Vn) |n)(n|, (12) 

Ti=nrnin Cp) 

nmaai(p)-l 

TVfe(-l)pTVTl(-l) =^p„+i,„+i sin^ (^Vn + l) \n){n\. 

Ti = max{0,n2nin (p) ~ 1} 

(13) 

By assumption, po € T7*. Then, O, ©-(HD above and 
induction on k show the assertions in Proposition [T] 


C. Computation of Qv{p,u) 

Fix any p G T7, and n G N. Recall that V (p) = 
Tr (d(JV)p), where cT; N —K be given by d(n) = (n — n)^. 
Note that ® implies that, for each u G {—1,0,1}, 

E [V{pk+i) I Pk= p,Uk = u] 

= Tr (d{N)Mg{u)pMl{ujj -fTr (d{N)Me(u)pMl{u^. 

(14) 












Take u = 0. From ([H])-® in Appendix iBl one has 

E [v{pk+i) I Pk = p,uk = 0] 

= Tr (d{N)Mg{0)pMli0)^ + Tr (d{N)Me{0)pMl{0)j 
= Tr (d(Ar) [Mg(0)pM^(0) + Me(0)pMl;(0)]) 

= TTid{N)p) = V{p). 


In particular, 

Qvip,0) = 0. (15) 

Now, take u = +1. Then, Gil above and Coli-GD in 
Appendix [B] provide that 

E [V{pk+i ) \ pk = p,uk = +1] 

= Tr (sin^ (fVAf + l)d(7V + l)p) 

+ Tr (cos^ (^VAl + l)d(Ar)p) . 

By summing and subtracting Tr (sin^ + l)d(Ar)p), 

E [V{pk+i ) I Pk = p,uk = +1] 


= Tr {d{N)p) 

+ Tr (sin^ [d{N + 1) - d{N)] p) 

= V{p) + Tr (sin^ [d{N + 1) - d(Ar)] p) . 

In particular. 


Qy(p,+1) = -Tr (sin2(f y]vTT)[d(iV + 1) - d(Ar)]p), 
= - Pnn [2(n - n) + 1] sin^ {^Vn+ l) . (16) 

nSN 


Finally, take u = —1. Using Gl above and G3- 
G1 in Appendix |B] E[U(pfc+i) | p^ = p^Uk = -1] = 
Tr (sin^^f ^)d(Ar - l)p) + Tr (cos2(^/V)d(iV)p). 

By summing and subtracting Tr ^sin^ d{N)p'^, 


E [u(pfc+i) \ Pk = p,uk = -1] 

= Tr {d(N)p) + Tr (sin^ Tiv) [d(Ar - 1) - d(iV)]p) 
= V{p) + Tr (sin^ [d{N - 1) - d{N)] p) . 

In particular, 

Qy(p, -1) = - Tr (sin^ (f Vtv) [d{N - 1) - d{N)] p) , 
= [2(n-n) - l]sin^ (^V^) . (17) 

nGN 


D. Proof of Lemma Q] 

Assume that Nq is even (otherwise one may take Nq + 1 
instead of Nq in this proof). Define the function p: N —?> M 
by 



By definition, one has + f for all £ G N. 

Let h = TV/A — arcsin ^•\/l/2 — . Using the definition of 


h and the symmetrie^ of the function sin^(-), it is easy to 
show that 


1/2 — a < sin^(a; + 7r/4) < a + 1/2, Va: £ [—/i, h]. (19) 


Let £ G N be even and big enough such that the following 
two conditions are simultaneously met: 

p(£) > iVo/2 + iV, \eQNQ/^p(l)-NQ/2<h. (20) 


Now, take N = [p(iy\ — "^ + 1 > N, where [rf] denotes the 
greatest integer which is less or equal to p. By construction, 
p(i) is in-between the points + Ao/2 — 1 and N + Nq/2, 
and hence it is in the interval [N, A^-fiVo — 1]. Then, for n = 
N,..., N+Nq — I, one has that |n—p(£)| < Nq/2. Consider 
the function (j){x) = ^y/x. From the fact that 4>'{x) = 
by the mean value theorem applied to the function f and the 
second inequality in (l20l i. one obtains 




for n = A^,..., A^ + A^o - 1- 


Then, the proof follows easily from Gl . Gl and the fact 
that sin^(a: — £7r/2) = sin^(a:), for every even £ G N. 


E. Proof of Lemma \2\ 

Proof of the first claim: Let u G {—1,0,1}, p G I?*. Recall 
that lU(p) = -EnGN^L- Since Tr (p) = EnGN^’^" = 
1, then -1 = ^ ^(/°) ^ 0- Now, 

by G), E[lU(p/c+i) \ pk = p,Uk = u] = Pg,fe 1U(pEi) + 
Pe,kW{pl^fj, where pg^k,Pe.k > 0 with pg^k + Pe,k = 1- 
Thus -1 < E[lU(pfc+i) I Pk = p,Uk = u] < 0. Since 
Qw{p,u) is the difference of two numbers that are in- 
between — 1 and 0, one concludes that \Qw{Ptu)\ < 1. 

The second, third and fourth claims, are immediate from 
Gl . GSJ and GtT i in Appendix O respectively. 


F. Proof of Lemma |5] 

Proof of the first claim: Let p G Dmg- By (l8]i-(l9]i in 

Appendix|B] Mg{Q)pM^g{Q) + Mf.{Q)pM\(ff) = p. Taking 
Pfe = p in Ufc = 0 in G). define 


Mg(0)pMt(0) 

Tr (Mg(O)pMi(O)) ’ 


for y = g,e. 


Hence, ap® -I- (1 — a)p® = p, where a = Pg^k = 
Tr (Mg(O)pMi(O)). In particular, ap9„ + (1 - a)pE = 

Pnn, for n G N. Note that, if a = 0, then Mg{0)pM^g{0) = 
0, and so p® = p. Similarly, a = 1 implies p® = p. Thus, 
the identity ap®„ + (1 — ct)pnn = Pnn, for n G N, still holds 
when q; = 0 or a = 1. From G), © and a = Pg^k, one has 


Qw{p,0) =W{p)- [pg,felU(pEi) + Pe,kW[pl+f)] 

= a {pin? + (1 - a) {pin? - [<^Pin + (1 “ o)pin] ^ 

riGN 

= a{l-a)Y,[pln-pln?>^, (21) 

nGN 


®More precisely, sin^(7r/2 — x) = 1 — cos^(7r/2 — a:) = 1 — sin^(a;). 


















thereby showing the first part of the first claim. 

If p = \m){m\ for some m S N with 0 < a < 1, then 
(|8]l-(|9]l in Appendix |B] imply that = p, and so 

Qw{p,0) = 0. Now, one shows that p = \m)irn\ for some 
TO £ N whenever Qry(p, 0) = 0. Suppose Qw{p,^) = 0. 
Then, (|2TI) implies that a = 0, or a = 1, or = p®^ for 
all n £ N with 0 < a < 1. Assume that a = 0. Hence, 
Mg{Q)pMl{Q) = EnGNP™cos2(^2^^)|n)(n| = 0 by 
® in Appendix [B] Suppose that p ^ |to)(to| for every 
TO £ N. Thus, there exists ni,n 2 £ N with ni ^ n 2 , 
Pm.m > 0, Pria.na > 0. Recall that sin(a;i) = ±sin(x 2 ) 
if and only if xi + X 2 = (-tt or X 2 — xi = (tt, where £ is an 
integer. Therefore, = ifc sin(which 

contradicts the assumptions that (j)o /tt is an irrational number 
and (pH = Tr/2 — npQ. One has shown that p = |to)(to| for 
some TO £ N whenever a = 0. If a = 1, or p®„ = p®„ for 
all n £ N with 0 < a < 1, then from similar arguments and 
computations one also concludes that p = |to)(to| for some 

TO £ N. 

Proof of the second claim: Let to £ N and take p = 
|to)(to| £ D*. It is clear that W{p) = = “I- 

From (fT0l i- (fT3T l in Appendix IbI one has that: 

(Mg(+l)pM^(+l))^^ = 6 {n, TO + 1) sin^ {^y/m+l) , 

{Me{+l)pMl{+l))_^^ = (5(n,TO)cos2 + , 

(Mg(-l)pMj(-l))^^ = S{n, to) cos2 , 

{Mei-l)pMl{-l)) = Sin + l,m)sm^ (tV™) , 

(22) 

where 6 {n,m) is the usual Kronecker delta: S{n,m) = 0 if 
n ^ m, and S{n,m) = 1 if n = to. In particular: 


Tr ^Mg(+l)pM^(+l)^ = sin^ {^Vm + l) , 
Tr (^Mei+l)pMl{+l)^ = cos^ (f , 

Tr (Mg(-l)pMt(-l)) = cos^ (f , 

Tr (Me(-l)pMt(-l)) = sin^ , 


E 

nGN 



My{u)pMl{u) 

(^My{u)pMliu)) 


2 


= 1, for u = ±1, y = g,e 

nn 


(assuming no division by 0). Now, using O and the above 
computations, one gets 


E [W (Pfc+l) I Pk= P:Uk= ±1] 

= Pg,kW{pl^^)+Pe,kW{pl+^) 

= - ^ [Tr(Mg(±l)pMt{±l)) X 

^ y f My{±l)pMli±l) \ 
nGN \Tr(Mg(±l)pMt(±l)) ) 
= -l = W{p). 


Therefore, Qvv(|w)(to|,±1) = 0. 


G. Proof of Proposition \2\ 

Fix p £ D,. Since Tr(d(Af)p) = Tr (d(7V)Ap) and 
Pnn = {£^p)nn foi' 71 £ N, the first two assertions are 
immediate from the definitions. As for the third and fourth 
assertions, let Ipj) = ^ 

p|to) = Pmn\n), for TO £ N. Using ©-(gli: 

Mg(0)pMt(0)|V>) 

ramax(p) 

m,n —0 

Me( 0 )pMt( 0 )|V-) 



m,n —0 


Myi + l)pMli + im 

nma.xip) + i 

= ^ prn-l.nSm{^^/m) sin{^VnTT) (m\'ip)\n + 1), 

m—l,n—0 

Me{+l)pMl{+l)\i;) 

rimxxip) 

= ^ Pmn cos {^Vm + l) cos (^Vn+ l) (TO|V')|n), 

m^n—0 

My{-l)pMl{-m) 

rimxxip) 

= ^ Pmn COS COS {m\tp)\n), 

m,n—0 

M,i-l)pMli-lM 

‘G.max (p) 

= X! Pm+i,nSin (^\/toTT) sin (^s/n) (TO|V')|n - 1). 

m—0,n—l 

Since Ap £ £), C D,, Umaxi^p) = rimaxip) and 
{Ap)nn = Pnn, the proof is straightforward from (l8Tl- (fT3T l 
in Appendix iBl 
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